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1. Introduction 

For over three decades lattice self-avoiding walks have been of interest both as models 
of polymers in dilute solution and as interesting and non-trivial problems in Statistical 
Mechanics [T]. The motivation for using these lattice models for the modelling of real 
polymers in solution comes from considerations of Universality; if the essential features 
are present in the minimal model, then it should accurately represent the critical 
behaviour of the real system. The essential features were identified as the excluded 
volume interaction and an effective attractive interaction modelling the difference in 
the solvent-monomer and monomer-monomer affinities. As the temperature (or solvent 
quality) is changed, the competition between these interactions gives rise to a collapse 
transition (the point) which separates the good solvent and bad solvent phases. 

Lattice walk models are coarse-grained representations of real polymers, and so the 
precise details of how these essential features are incorporated should not matter. Whilst 
the standard interacting self-avoiding walk (ISAW) model, where walks are forbidden 
from visiting a lattice site or lattice bond more than once, is the canonical model to 
study polymers in dilute solution, two other models were presented as alternatives: the 
vertex- interacting self- avoiding walk (VIS AW), and O(n=0) symmetric walk introduced 
by Blote and Nienhuisp] where the walk is allowed to visit sites twice, but not cross 
itself, and the interacting self-avoiding trail (ISAT), where the walk is allowed to visit 
sites twice and cross[5]. The self-attraction is included between non-consecutive nearest- 
neighbour visited sites for the ISAW, but is associated with the doubly visited sites in 
the other two models. 

Simple universality arguments would lead one to think that these models should be 
in the same universality class, both in good solvent and at the collapse transition. Whilst 
this seems to be the case in good solvent, exact results for the two-dimensional ISAW 
and the VISAW models show that these two models are not in the same universality class 
at the collapse transition, the first having a value for the thermal exponent vg = 4/7 [3] 
whilst the second has vg = 12/23 [4|. 

The situation for the 2D ISAT is far less clear; for the moment there are no exact 
results, but a wide range of estimates for vg. In the eighties the ISAT at the collapse 
point was in a different universality class than the ISAW[6], whilst in the early nineties 
some authors claimed to find evidence that the two were in the same universality class[7]. 
In 1995 Owczarek and Prellberg|8j studied a kinetically growing self-avoiding trail model 
with no interaction. This model may be mapped onto the ISAT with a particular value 
of the attractive interaction. They found a value of v = 1/2. This result could lead one 
to conclude that the kinetic self-avoiding trail maps onto the ISAT in the bad-solvent 
regime. They exclude this possibility by showing that the density of the walk vanishes in 
the infinite walk limit. In 2007 Owczarek and Prellberg[9] confirm some of their results 
with a direct FlatPERM simulation directly on the ISAT model for walks up to about 
2 000 000 steps. 

In this article we re-examine the ISAT model using a numerically exact transfer- 
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matrix calculation in the full fugacity/inter action plane. We give compelling evidence 
that, contrary to previous claims, the ISAT model in two dimensions is in the same 
universality class as the VISAW, with a thermal exponent ug = 12/23. This is reinforced 
by the presence of a phase transition line separating two finite-density phases which 
we conjecture to be in the Ising universality class, also present in the VISAW phase 
diagram [2]. 

This paper is organised as follows: the ISAT model is presented, followed by 
the results obtained from the transfer matrix calculation. The article ends with a 
discussion of possible reasons for the apparent difference of results between those found 
by Owczarek and Prellbergpj], and those found here, and their consequence for the study 
of self- avoiding walk models where frustration effects become important. Such models 
are of increasing interest as toy models for biopolymers [1 Oj . and as such it is important 
to understand in detail the effect the underlying lattice has on the critical behaviour of 
the model, and under what conditions such a competition may arise. 



2. Model and transfer-matrix method 



The ISAT model studied here is defined as follows: consider all random walks on the 
square lattice which do not visit any lattice bond more than once. Doubly visited sites 
may correspond to either crossings or "collisions" , both are assigned an attractive energy 
—e. The partition function for the model is 

z=Y. rNtNi > (!) 

walks 

where K is the step fugacity, r = exp(/3e), N is the length of the walk, and N I is the 
number of doubly-visited sites. 

This partition function may be calculated exactly on a strip of length L x — > oo and 
of finite width L by defining a transfer matrix T. If periodic boundary conditions are 
assumed in both directions, the partition function for the strip is given by: 

Z L = lim Tr (T l A . (2) 

The free energy per lattice site, the density, and correlation length for the infinite strip 
may be calculated from the eigenvalues of the transfer matrix: 

/ =^ln(A ), (3) 

(XT \ K dX ° (A\ 




e(X,r)= (ln^ ) , (5) 

where Ao and Ai are the largest and second largest (in modulus) eigenvalues. 

It is expected that Z, p and £ should have the following scaling forms close to the 
critical fugacity (for fixed r): 

Z ~\K-K c \-i, (6) 
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£ ~\K-K e \- V , (7) 
p L (K) = Poo (K) + LV»- 2 p(\K - K c \L l ' v ). (8) 

Z corresponds to the high temperature expansion of the susceptibility of an equivalent 
magnetic model, hence the use of the exponent 7. 

These scaling properties enable estimates of the critical lines to be calculated using a 
phenomenological renormalisation group method. For example a critical point estimate 
for a pair of lattice widths L and L' is given by the solution of the equation: 

T = l < 9 > 

with estimates of the critical exponent v given by: 

1 _ 1U & \dK I dK J 1 nm 
log (L/L') • 1 1 

The critical dimensions of the magnetisation and energy fields may be calculated from 
the first few eigenvalues of the transfer matrix: 

Lin \ 



1 n 



2vr 
Lin 
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(12) 



2tt 1 

The scaling dimensions x a and x e may be related to the correlation length exponent u 
and the exponent 7 through standard relations 

^=^— , (13) 
7 = 2z/(l-x CT ). (14) 

For a more detailed discussion of the transfer matrix method, the reader is referred 
to the article of Blote and Nienhuis [2]. 



3. Results 



The transfer matrix for a lattice walk breaks down naturally into three sectors: the 
empty lattice sector (a 1 by 1 block), and two sectors corresponding to an even or odd 
number of horizontal links on a lattice column. In the zero-density phase, the largest 
eigenvalue is given by Ao = 1, corresponding to an empty lattice. In the dense phase one 
may take the largest and second largest eigenvalues from different sectors. For walks on 
an odd lattice width, the largest eigenvalue, A D , of the odd sector is always larger than 
the largest eigenvalue, A e , of the even sector. For even lattice sizes there is a line in the 
(K, t) plane where A Q = A e . A crossing of the two largest eigenvalues indicates a critical 
line. Such a crossing is not normally expected for a finite lattice width, but occurs in such 
walk models, and often indicates a transition between a crystalline phase and a liquid 
phase. The existence of such a phase transition is corroborated by phenomenological 
RG. The phase diagram calculated for even lattice sizes is shown in Figure [H 
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Figure 1. Phase Diagram calculated using Phenomenological RG (equation [9|) for 
even lattice sizes using matrices in the sector with the largest eigenvalue. The upper 
line is estimated with the crossings of the eigenvalues between even and odd sectors, 
(colour online) 



The phase diagram shows three phases: the zero-density phase, a crystalline phase 
and a liquid phase. This phase diagram is different from the phase diagram for the 
ISAW model for the point, where there is only one high-density phase. The phase 
diagram is qualitatively similar to the phase diagram of different models which display 
frustration effects due to a competition with the underlying square lattice. In such 
models the details of the critical behaviour on the crystal/liquid phase transition and of 
the multi-critical point at coexistence between the three phases depend sensitively on 
the details of the attractive interaction 

An estimate of the location of the multicritical point may be found from the 
crossings of the estimates of v as a function of r, shown in figure [2] for odd lattice 
widths. In the infinite lattice limit v = 3/4 for r < Tg and v = 1/2 for r > Tg, adopting 
a non-trivial value for r = Tg. The lines cross at or very close to vg = 12/23 ~ 0.52174, 
which is the exact value for the equivalent point in the VISAW[4], and far from previous 
conjectures of vg = 4/7 ~ 0.57143[7J (universality class of the ISAW model) or v = 1/2 
proposed by Owzcarek and Prellberg[8J. What is interesting for odd lattice sizes is 
that for r = 3 (the conjectured location of the collapse transition in this model) 
all the solutions of the phenomenological RG equation [9] occur at K = 1/3 with 
A = Ao = 1. A , the largest eigenvalue of the odd sector, corresponding to the second 
largest eigenvalue of the transfer matrix, is smaller than 1 for both r < 3 and r > 3. 
This singular behaviour means that the derivative needed in Equation ( flOl) is undefined, 
and the estimate for v exactly at r = 3 is missing. The results for even lattice sizes are 
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Figure 2. Estimates of v from equations [5] and [TO] for odd lattice sizes. The horizontal 
line corresponds to v = 12/23. In the intersection region there is a data point every 
0.01 along the x-axis, except at the point t = 3.00 (see text), (colour online) 

given in table [TJ These results are consistent with those found using odd lattice sizes, 
as well as an alternative phenomenological RG based on the scaling of the density. The 
value of Kg and Tg converge nicely to the values Kg = 1/3 and Tg = 3. The estimates of 
vq, whilst remaining close to the expected value of 12/23, they overshoot. It sometimes 
occurs that estimates overshoot their asymptotic values, reaching a maximum before 
converging, and is already the case for the SAWjl]. With the limited number of lattice 
widths available here we do not see a maximum. To try and confirm this possibility, a 
different way of estimating vg is used. There are strong reasons to believe that Tg = 3 
corresponds to the collapse transition [5]. If at this point we find a value of v different 
from 3/4 and 1/2, this point must then be identified with the collapse transition, this 
was also the argument used in reference [5]. We calculate x e at fixed r = 3 using 
Equation (fl2l) with K solution of Equation ([9]). This gives us vg = 1/(2 — x £ ). Since 
two lattice widths are required to calculate K c {r = 3), this gives two estimates for ug, 
which are shown, along with estimates of Kg and x a , in table [2j These estimates of vg 
also overshoot 12/23, but they reach a maximum and seem to converge to the expected 
value. The small number of lattice sizes does not permit a fuller finite-size scaling 
analysis, but the different results presented seem to clearly support the identification 
of vg = 12/23, corresponding to vg for the VISAW, for which the exponent has been 
determined exactly [4]. We confirm the previous conjecture that the collapse is likely 
to occur at exactly Tg = 3 [8]. The numerical results for Kg are consistent with the 
identification Kg = 1/3. 
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Table 1. Results for the multicritical values of K, r and v calculated for even lattice 
sizes using phenomenological RG (equations 191 and fTU|) . The last line conjectures exact 
values for these parameters. The value given for v corresponds to the exact value for 
the VISAW, the value of r is the value of r for which the model maps onto a kinetically 
growing SAT, and the value of K is conjectured from the numerical results given here. 
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Table 2. Results for Kg, x a and v% calculated using phenomenological renormalisation 
group, fixing tq — 3. The thermal exponent is calculated via the scaling exponent x e 
and using Equation (fl"3)) . For each point calculated, there are two values of x e , one for 
the smaller lattice width L and one for the larger lattice width V . 



L/L' 


K e 


X a 


u e = 1/(2 -x s (L)) 


u e = 1/(2 - x e (L')) 


2/4 


0.333865 


0.078111 


0.520577 


0.521291 


4/6 


0.333259 


0.085770 


0.522540 


0.522700 


6/8 


0.333221 


0.086325 


0.523004 


0.523335 


8/10 


0.333246 


0.085686 


0.522984 


0.523307 


10/12 


0.333269 


0.084817 


0.522830 


0.523118 


conjecture 


1/3 


1/12 = 0.08333- •■ 


12/23 = 0.521739- ■• 


12/23 



The VISAW model also displays a liquid/crystal phase transition, found to be in 
the Ising universality class[2]. If the collapse transition is of the same type here as for 
the VISAW model, the liquid/crystal phase transition here should also be in the Ising 
universality class. The exponent values have been calculated for odd and even lattice 
sizes. Due to parity effects, the odd and even lattice sizes give two lines of estimates, 
both of which converge (one from above, the other from below) leading to v = 1.00±0.03, 
consistent with an Ising universality class. 

All the thermal exponents seem to coincide with those for the VISAW model. We 
also calculated the magnetic critical dimension x a ~ 0.083 ± 0.002 (to compare with 
1/12 = 0.0833333). If x a = 1/12 and v = 12/23, then le = 22/23. This is different 
from the VISAW model for which x a = —5/48 (or 70 = 53/46) [4J. This difference 
reflects the larger configuration space opened up by allowing the walk to cross at sites. 
Similar differences are seen between the ISAW model on the square lattice and on the 
Manhattan lattice [T2]. 

The density at the collapse transition for the ISAT is shown in Figure [31 At first 
sight it seems to indicate a finite density for the infinite system, but when it is fitted 
with the scaling relation[8l an excellent fit is found for = if we use v = 12/23. We 
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Figure 3. Density calculated at r = 3 setting Ai(i^) = 1 for odd lattice sizes. The 
solid line represents a fit to the scaling law p(L) = p 00 J r aL 1 l v ^ 2 with p^ = 0, a = 1.65 
and v = 12/23. (colour online) 

were not able to fit with ^=l/2orz/ = 4/7, however, given the number of data points, 
and small lattice widths examined, it cannot be excluded that other good fits could be 
found for other the exponent values if additional correction terms are included. It is, 
however, a reassuring consistency check, and indicates that our results are consistent 
with the claim of Owczarek and Prellberg[8] that the density is indeed zero in the infinite 
walk limit. 

4. Discussion 

In this paper results indicating that the ISAT model at the collapse transition is in the 
same class of universality as the VIS AW model introduced by Blote and Nienhuis[2j are 
presented. The correlation length exponent is consistent with vg = 12/23. These results 
are at variance with previous results, most notably of Meirovitch and coworkers [7 J who 
conjectured that the model was in the same class as the standard ISAW model, and 
Owczarek and Prellberg who give the correlation length exponent as vg = 1/2[8J. In the 
first case, the model was studied using the scanning Monte-Carlo method. It is known 
that the calculated critical exponents are sensitive to the estimations of the location 
of the multi-critical point, and their estimated critical point, whilst close to ours, is 
significantly lower (r B = 2.962 ± 0.004)0. 

The apparent contradiction with the results of Owczarek and PrellbergjH] is more 
interesting. They performed Monte-Carlo simulations for extremely long chains at the 
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same value of tq = 3, and claimed to find v = 1/2, clearly excluded from our results. 
However, they used the often-used identification of the exponent v with the radius of 
gyration: 

(Rg) ~N v . (15) 

This equation defines v as a geometric exponent, equal to the inverse of the Hausdorff 
fractal dimension of the walk. When the polymer is collapsed, v = 1/2 (in two 
dimensions). This occurs along the first-order line separating the zero-density phase 
and the crystalline phase, but the thermal exponent v is not defined here, since there is 
no diverging correlation length. That equation[15]is not always valid is trivially apparent 
along the liquid/crystalline transition, where the dimension of the walk is 2, but the 
exponent v — 1. We suggest that Owczarek and Prellberg have correctly identified the 
Hausdorff dimension of the walk to be du = 2, but that once the dimension of the walk 
and the lattice are the same, equation [15] no longer applies. Since the polymer is "space 
filling" (even if in this particular case = 0), it "sees" the underlying lattice, allowing 
for competition between the short-range interactions and the lattice geometry. We 
believe this to be the origin of the difference between the ISAW and both the ISAT and 
VISAW models[13], and the apparent lack of universality in these lattice walk models. 

The connection between the VISAW model and the ISAT model needs to be further 
investigated, and the particularly nice values of Kg = 1/3 and r e = 3 leads one to ask 
if an exact resolution of the problem would not be possible. In any case, as for any 
numerical calculation, an independent verification of these results by other methods 
would be welcome. 
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